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Recently, several interesting proposals were made modifying the law of gravity on large scales, 
within a sensible relativistic formulation. This allows a precise formulation of the idea that such a 
modification might account for galaxy rotation curves, instead of the usual interpretation of these 
curves as evidence for dark matter. We here summarize several observational constraints which any 
such modification must satisfy, and which we believe make more chaUenging any interpretation of 
galaxy rotation curves in terms of new gravitational physics. 



PACS numbers: 

I. INTRODUCTION 

The observed flatness of galaxy rotation curves are 
used as important evidence for the existence of dark mat- 
ter on galactic scales. These observations thus form a 
keystone of the structure of inferences that produce the 
current cosmological paradigm, which requires the vast 
majority of the matter in the universe to consist of un- 
known, dark forms. Although this paradigm has received 
considerable support from recent observations, this sup- 
port has been at the expense of simplicity, since the ev- 
idence now requires two forms of dark matter: both the 
conventional sort plus the newly-discovered 'dark energy' 

m- 

The validity of the Newton-Einstein laws of gravity on 
galactic scales is an important component of the stan- 
dard reasoning, and one for which the direct observa- 
tional support is comparatively weak. Although much 
less work has been done to explore possible modifications 
of gravity on large scales than has been invested in dark 
matter models, over the years a literature on the subject 
has emerged Q. Among the variants is a reason- 
ably simple modification of Newtonian dynamics which 
describes galaxy rotation curves quite well, without re- 
quiring the existence of dark matter 11-1^]. 

Much of this work has remained outside of the main- 
stream for various reasons. Some proposals invent ad- hoc 
nonrelativistic potentials without addressing how these 
potentials might arise within the context of relativistic 
field theory. This is not a small omission, because it is 
precisely the issue of relativistic consistency which makes 
modifying physics at large distances difficult. Although 
it has long been known that the laws of gravity are likely 
to receive modifications on short distance scales, general 
constraints like unitarity and lorentz invariance in weak 
fields limit what can be done on very large scales. These 
arguments seem to restrict any such modification to take 
the form of supplementing the graviton with fields de- 
scribing various low-spin degrees of freedom (l5| , leading 
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to scalar-tensor-Maxwell kinds of theories at long dis- 
tances.]^ While some of the proposed modifications to 
gravity fall into this category it has proven difficult to 
reproduce the galactic rotation curves with these theo- 
ries without also running afoul of other constraints |p^ . 

The situation may now be changing, with several re- 
cent claims that self-consistent particle-physics models 
permit the law of gravity to be modified on large scales. 
Several recent proposals do so using constructions 0] 
within the braneworld picture [ p^|Jl9|] , although none of 
these has yet tried to achieve acceptable phenomenology 
at large distances. Ref. aims to describe spiral galaxy 
rotation curves assuming the existence of extremely light 
four-dimensional fields trapped onto galactic-sized de- 
fects. These new developments motivate re-examining 
the astrophysical constraints on modifications to gravity. 

It is not our purpose in this note to present another 
modification of gravity at large distance scales. Rather, 
we put aside general issues of principle, and indicate what 
targets proposed models must hit in order to agree with 
the extensive astrophysical data on galaxy dynamics. We 
shall argue here that the data pose a significant challenge 
which must be met by any hypothetical modifications 
to gravity, and we believe none of the extant proposals 
completely meet this challenge. 

The constraints come in several types. First, the sys- 
tematics of galaxy rotation curves are observationally 
found to be described by a universal shape which de- 
pends only on the galaxy luminosity and visible size. In 
particular, the huge range in sizes of observed galaxies 
precludes theories which modify gravity at a fixed dis- 
tance for all galaxies. Second, the galactic gravitational 
potential for spiral galaxies like the Milky Way is mea- 



*These do not include higher-derivative theories, which 
generically have runaway modes which violate the unitarity 
and stability constraints. 
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sured in 3 dimensions, not just in the plane of the galactic 
disk. Finally, additional constraints come from the TuUy- 
Fisher relation and the x-ray observations of gas in galaxy 
clusters. Some of these constraints have been raised else- 
where as objections to various specific proposed modifi- 
cations 

We emphasize that our purpose in making this argu- 
ment is not to argue that phenomenologically-successful 
modifications to gravity cannot be made (although this 
may in the end prove to be true). Rather, we intend to 
provide a convenient summary of the astrophysical con- 
straints which are relatively unfamiliar to the particle 
physics and quantum gravity communities. 



II. CONSTRAINTS 

We now summarize the most serious astrophysical con- 
straints, whose successful description we claim to be a 
necessary (though not sufficient) condition for the suc- 
cess of any putative realistic modification to gravity at 
large distances. Since most of the proposed modifications 
are motivated to describe galaxy rotation curves, we de- 
scribe these first before turning to the other constraints. 

1. Rotation Curve Systematics: Observations on many 
galaxies indicate rotation velocities, v^ot, which for suf- 
ficiently large distance, r, from the galactic center, be- 
come independent of r. This is to be contrasted with 
the Keplerian form, u^^j ~ GM/r, which would be ex- 
pected for circular orbits well outside of a concentrated 
mass distribution, and which Newton's Law would also 
approximately predict for galaxies sufficiently far outside 
the visible material (in the absence of dark matter). 

A common way to account for this amongst proposed 
modifications of gravity is to have the gravitational po- 
tential be modified from the Newtonian form —GM/r to 
a new form: —GM/r + U{r). (Proposals for U{r) in- 
clude linear logarithmic |§,go| or Yukawa-type 0] 
functions of r, or have it interpolate between different 
asymptotic Newtonian regimes with different values for 
G The successful comparison of General Relativ- 
ity with solar-system and terrestrial measurements is 
then accounted for by ensuring U(r) becomes sufficiently 
small when r < ro, where tq is a new fundamental char- 
acteristic distance scale. 

For the vast majority of models, the modified force 
is itself proportional to AI - i.e. U{r) = M u{r) with 
u{r) an Af- independent function - because it is obtained 
by the superposition of forces exerted by the individual 
particles of which the galaxy is made. It follows that 
any such model crosses over from Newtonian to exotic 
force laws for a radius, ro, which is independent of the 



properties of the galaxy involved. ^ 

Any model for which this crossover radius is universal 
- regardless of the details of the force law responsible - 
does not provide a good description of the systematics of 
galaxy rotation curves. It does not because it predicts 
the scale, rg, to be independent of the luminosity of the 
galaxy. This conflicts with observations, for which the 
transition radius (at which the rotation curves deviate 
from the predictions of Newtonian gravity applied to the 
visible mass) occurs at radii which are correlated with 
the galaxy's optical size,^ -Ropt- Since 'non-Newtonian' 
rotation curves have been seen in galaxies that range over 
several orders of magnitude in overall luminosity, the ob- 
served transition radius (and optical size) varies over an 
order of magnitude, from a few to a few tens of kilopar- 
secs. 

The data indicate that when a galaxy's rotation speed 
is plotted against r/i?opt, the resulting curve is a univer- 
sal function which depends only on the total luminosity 
of the galaxy (see Eq. (14) and Fig. 10 in (HJ). In galax- 
ies with low luminosities, the contribution of dark matter 
becomes important well inside the optical radius, while 
in highest luminosity (very large) galaxies luminous mat- 
ter adequately describes the rotation curves to at least 
r/Ropt ''^ 1. This behavior is incompatible with modifi- 
cations of the gravitational force law at a fixed distance 
from the source. 

To be completely concrete, consider the case of two 
galaxies, N598 and N801 (with data taken from p3[). 
The data show that the rotation curve of the low- 
luminosity galaxy, N598, would demand the modification 
of gravity at distances less than its optical radius, which 
is 3.9 kpc. This would cause a dramatic effect on the 
rotation curve of the high- luminosity N801, well inside 
of its optical radius of 38.4 kpc. No such distortion is 
observed. We see that the larger galaxies would exhibit 
more - not less ~ dark matter than the small ones at 
a given fraction of their optical radius, were the law of 
gravity modified at a fixed distance. 

The existence of modifications to gravity which can 
successfully ensure that the modifications become impor- 
tant at distances which depend on the galaxy involved 
is demonstrated by the proposal called MOND [p]-p^. 
MOND's success in describing galaxy rotation curves re- 
lies on its modifying Newtonian dynamics at a minimum 
acceleration, rather than modifying Newton's Force Law 
above a maximum distance. That is, the modification 
of rotation curves from the Newtonian prediction does 



^Ref. avoids this conclusion by adding two additional 
force terms, one depending on M and one depending solely 
on the distance from the galaxy center. 

*We use i?opt = 3.2Rd, where Rd is the exponential disk 
scale length, in accordance with the definition of Ref. [|23| 
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not happen at the same distance from the center of all 
galaxies, but occurs at a radius for which the Keplerian 
acceleration is of order ag ~ 10^^" m/s'^ ^4 24|. Indeed, 
part of what is remarkable about the observational suc- 
cess that MOND has in describing observations is the 
successful fitting of the rotation curves of many different 
kinds of galaxies in terms of a single function describing 
the crossover from Newton's law at high acceleration a 
to its modified version at small acceleration around oq. 
Best fits also give very reasonable values for the mass-to- 
luminosity ratio, T = M/L (and so does not introduce a 
new reason to require a new kind of dark matter). 

2. The Tully-Fisher Relation and M/L: Besides pre- 
dicting flat rotation curves, any modification of gravity 
should also account for the dependence of the rotation 
curves on galaxy luminosity. Part of this dependence is 
summarized by the Tully-Fisher relation which ex- 
presses an observational fact: the luminosity, L, of a spi- 
ral galaxy is tightly correlated with the asymptotic limit 
of its rotation velocity, Voo = lim,,^oo Wrot(?'), according 
to 



L oc v^, 



(1) 



with a « 4 |2|]. 

While the Tully-Fisher relation (and particularly its 
small scatter) is not completely understood within the 
dark-matter picture, it can be roughly derived in sim- 
ple models where galactic disks collapse with conserved 
angular momentum within dark matter halos of a form 
predicted by numerical simulations p7[ |. But any pro- 
posal to replace dark matter with a modification of grav- 
ity should be able to do better, by directly yielding the 
T-F relation once the mass-to-light ratio T of the matter 
constituting the galaxies is specified. Conversely, to sat- 
isfy the Tully-Fisher relation any modified gravity theory 
will predict a specific form for T(L), and this prediction 
can be tested. 

For example, all models which modify the position- 
dependence of the gravitational potential at large r by 
supplementing the Newtonian term with M u(r) predict 
v^/M to be purely a function of r for circular orbits. If 
u{r) for a given model is designed to approach a con- 
stant for large r (to describe the observed flat rota- 
tion curves), then that model predicts M oc p2| , 
which is consistent with the Tully-Fisher relation only if 
T — M/L (X L~^/^. This relation is, however, strongly 
constrained because T can be estimated both theoreti- 
cally and observationally for baryonic matter in galaxies 
over a wide range in luminosity. 

Theoretically, T « (0.5 — A)Mq/Lq in stellar ma- 
terial is predicted by stellar evolution models that can 
reasonably reproduce the colors and spectra of observed 
spiral galaxies Similar values can also be directly 

measured in the centers of both high- and moderate- 



luminosity galaxies, where stars are expected to domi- 
nate the mass (and where Newtonian gravity would pre- 
sumably hold) . More generally, stars are observed to 
dominate the mass of known baryonic matter in galaxies 
- only in the smallest dwarf galaxies does the mass in gas 
approach or (marginally) exceed that in stars | |30[ | . 

Thus T ^ (0.5— 4:)Mq/ Lq is expected to hold for spiral 
disks of all luminosities. This does not pose a problem for 
MOND (for example) , which predicts L oc and there- 
fore generally predicts T to be independent of L and 
to lie within this range |3H |. It is problematic for any 
theory predicting a varying T(L), since the T-F relation 
is observed to hold over four decades in L For ex- 
ample, the prediction T oc L~^/'^ generic to many models 
implies that T must be 100 times larger in the least lumi- 
nous observed galaxies than in the most luminous ones. 
Since the stellar populations of dwarf and giant galax- 
ies probably have similar T, the T-F relation can be fit 
only by assuming that ~ 99% of the mass of the least 
luminous galaxies is in an unknown form (i.e. neither 
observed stars or atomic or molecular gas). Although it 
may be possible to invent forms of baryonic matter in 
which such large amounts of mass might be hidden, this 
would be rather unsatisfying in a model constructed to 
remove the need for dark matter. 

3. Galaxy Potentials in Three Dimensions: Since New- 
ton's law in d space dimensions implies a gravitational 
potential which varies as (l/r)''^^, it is tempting to try to 
obtain logarithmic potentials - and so constant Keplerian 
velocities - by trapping the field which mediates this force 
in some way to two space dimensions (for which Newton 
predicts a logarithmic dependence). At first sight this is 
a very attractive possibility, since if a field contributing 
to long-distance interactions is trapped on a surface of 
thickness £ ~ several kiloparsecs, then interparticle po- 
tentials vary like l/r when r <^ £, but cross over (within 
the plane) to logr dependence for r ^ £. (See ref. l2C| ] 
for a recent attempt to construct such a mechanism.) 

A fundamental prediction of any such mechanism is 
that galactic rotation curves should look flat in the direc- 
tions along the trapping surface, but the galactic grav- 
itational well should be Newtonian when examined in 
directions perpendicular to this surface. Putting aside 
questions of how elliptical galaxies (which are not flat- 
tened and yet also appear to have dark matter [ ^ ) 
might fit into such a scenario, there are several pieces 
of evidence concerning our own galaxy which prospective 
models must confront. 

• Globular Clusters: Globular clusters are well-defined 
clusters of stars which are distributed in an approxi- 
mately spherically-symmetric way about the galactic cen- 
ter. They provide an important constraint on models 
of the galaxy because: (i) they are not confined to the 
plane of the spiral arms; (m) they are distributed out to 
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40 kpc from the galactic center; and (iii) at least 26 of 
the roughly 150 known globular clusters have measurable 
proper motions relative to extragalactic background ob- 
jects [ p3|j3^ ], and so have known velocities and positions. 

The constraint on galaxy models arises as follows. 
Given their known positions and velocities, it is possi- 
ble to reconstruct the galactic orbits of many globular 
clusters given a model of the galactic gravitational po- 
tential. In particular, for many models - including those 
for which the galactic mass does not extend out to 40 
kpc or so - the orbital apocenters {i.e. those points most 
distant from the galactic center) lie 50 to 80 kpc from 
the galactic centers. Since most of the time in an orbit 
is spent near the orbital apocenter, it is statistically very 
unlikely to find less than a handful of clusters further 
than 40 kpc from the galactic center, as is observed 

• The Magellanic Stream: The Magellanic Stream is a 
trail of neutral hydrogen which extends in a great circle 
more than 100° across the sky, starting at the Magellanic 
Clouds. This hydrogen is most likely gas which has been 
stripped from the Clouds due to tidal forces, and because 
the Milky Way's gravitational potential is much deeper 
than that of the Clouds, the Stream is distributed along 
the orbit of the Magellanic Clouds themselves. This pic- 
ture is consistent with the measured infall velocity of the 
Stream, which varies nearly linearly as a function of dis- 
tance from the Clouds, reaching —200 km/s at the end 
furthest from the Clouds. The line-of-sight velocity of 
the center of mass of the Magellanic Clouds themselves 
is 61 km/s, indicating that the Clouds are probably near 
the pericenter - closest approach - of their orbit, which 
must be quite eccentric. 

The constraint on the Galactic potential is obtained 
by requiring the existence of an orbit with a pericenter 
of roughly q = 50 kpc and on which radial velocities as 
large as 200 km/ s are obtained |Q . This puts strong con- 
straints on models with all of the Milky Way's mass lo- 
calised where the luminous matter is, and in which grav- 
ity is Newtonian out of the galactic plane, since for these 
the Stream motion should be approximately Keplerian. 
Since the maximum radial velocity of a bound Keplerian 
orbit with pericenter q is Wmax — {GM/2q)^ , for q — 5Q 
kpc, the observed 200 km/s infall of Stream material im- 
plies a galactic mass A/ > 9 x 10^^ solar masses. This is 
an order of magnitude larger than the mass obtained by 
counting the luminous matter in the disk, but is consis- 
tent with the amount of dark matter required to account 
for rotation curves. 

• Local Group Galaxies: Constraints similar in spirit to, 
but weaker than, the one just described can be obtained 
from the dwarf galaxies which make up much of our local 
galactic neighborhood. A number of arguments indicate 
that these are bound to our galaxy and their observed 
distances and radial velocities constrain the shape and 
depth of the galaxy's gravitational well. In particular, the 



large speed of some of these dwarfs (such as Leo I, Pal 14 
and Eridanus) can only be bound by the presence of much 
more mass than is visible if the out-of-plane gravitational 
field is Newtonian. 

4-. The Dynamics and Structure of Clusters of Galaxies: 
As first notice by Zwicky [Q, the visible mass in clus- 
ters of galaxies is entirely too small to bind them, given 
the observed velocity dispersion (several hundred km/s) 
in their component galaxies; this was the first strong evi- 
dence for dark matter (or modified gravity) in cosmology, 
and is yet stronger today. X-ray measurements of clusters 
have revealed that their observable baryonic masses are 
dominated by hot gas of temperature ~ 1 — 10 keV, which 
also implies (using standard gravity) a binding mass well 
above that observed in galaxies and the hot gas itself. 

It is instructive to derive this inconsistency explicitly, 
in order to see the logic of the constraints posed on dark 
matter and on modifications of gravity. Consider, then, 
the equation of hydrostatic equilibrium in a spherical sys- 
tem (which describes intracluster gas well): 



IdP 

p dr 



(2) 



where p and P are the density and pressure of the gas, 
and a is the inward gravitational acceleration at radius 
r. Using the ideal gas equation of state, this can be 
re-written in terms of the gas temperature T and mean 
molecular weight, fi 0.6 (in units of the proton mass, 
TOp): 



dlogp dlogT _ r / pmp \ 
dlogr dlogr ^ T \ k J 



(3) 



where k is Boltzmann's constant. Clusters are observed 
to have temperature profiles, T{r), which are roughly 
constant outside of their cores. The density profile of the 
observed gas at large radii roughly follows a power law, 
p oc r", with index — 2 < a < —1.5 [Q. Using these, the 
fact that the gas dominates the cluster's baryonic mass, 
and Newton's expression for a(r) in the absence of dark 
matter gives 



fcr« (1.3- 1.8)keV 



Mr 



IMpc 



(4) 



for the baryonic mass within the central 1 Mpc of a typ- 
ical rich cluster, where M^ is the mass enclosed within 
r. The disparity between this estimate and the corre- 
sponding observed temperature ~ lOkeV indicates 
the need for dark matter or modified gravity. 

Similar reasoning provides a strong constraint on any 
modifications of gravity, for which the only matter 
present is visible and a{r) is calculable. Most generally, 
Eq. (^) implies an approximate relation between the to- 
tal mass, outer radius, and average temperature of the 
cluster. For a given gravity theory, this relation can be 
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checked by comparing to observed masses, temperatures, 
and sizes of observed clusters, as we have done for New- 
tonian gravity above. A useful compilation of applica- 
ble cluster data is given in |3|]. MOND can apparently 
survive this comparison as can the standard CDM 
cosmology with 5 — 10 times as much dark matter as 
baryonic matter, whereas (for example), the construction 
of would not account for the high cluster tempera- 
tures (and hence would still require the standard amount 
of dark matter) because most cluster gas would not lie 
within the domain walls and would hence feel only the 
usual Newtonian gravitation of the visible matter. 

A stronger constraint can be deduced using observed 
density, mass and temperature profiles of individual clus- 
ters. If gravity depends only on the visible mass, the 
observed density profile p{r) and enclosed mass profile 
Mr of a given cluster can be used to directly predict its 
temperature profile T{r) in a given gravity theory using 
Eq. (|^), and the profile compared to the measured one. 
Roughly speaking, the combination ra(r), which is pro- 
portional to T(r) (since din p/dlnr and T are observed 
to be approximately constant) should be nearly indepen- 
dent of r, as T(r) is observed to be. This technique is 
illustrated in detail in , where it is shown that MOND 
predicts rising radial temperature profiles that disagree 
badly with the observations (which show that clusters are 
roughly isothermal). 

Other modifications of gravity may encounter similar 
difficulties. Consider, for example, the hydrostatic equi- 
librium of cluster gas when large-r gravity is dominated 
by a linear potential PJlC|]. Since linear potentials are 
sensitive to the mass distribution at large radius, we will 
assume the cluster gas to not extend beyond an outer 
radius, For simplicity of analysis we bracket the 

real mass density for r < R between two extremes: it 
is taken to be concentrated at radius r = 0, or to be 
uniform within R. Under these assumptions we find that 
for r <^ R, ra{r) respectively grows either linearly or 
quadratically with r. Provided only that a more real- 
istic treatment lies between these two extremes, we see 
that T(r) is predicted to increase linearly to quadrati- 
cally with r. Such a prediction would be ruled out by 
the observations, although a more detailed study would 
be required to make this conclusion fully quantitative. 

Finally, we note that cluster mass estimates using 
galactic dynamics or the X-ray gas are consistent (to 
within a factor of two) with masses determined by weak 
gravitational lensing of background galaxies around the 
clusters 40 1. Thus any alternative gravity theory must 



^ Given the extreme sensitivity of the interaction energy to 
matter at arbitrarily large radii, it is unclear how to compute 
the characteristics of a given physical system with a linear 
potential without neglecting this contribution. 



also predict gravitational lensing in a way that allows 
this. This constraint was used, for example, by Q to 
rule out their model, and provides difficulties for MOND 
in cluster cores [|4|. 

It would be wonderful to discover new gravitational 
physics from astrophysical observations, and it is intrigu- 
ing that current braneworld models seem to imply that 
long-distance gravitational physics can be much richer 
than had been previously thought. In order for any such 
physics to provide a viable alternative to dark matter it 
must satisiy the above constraints, and a theory that suc- 
cessfully does so would be of sufficient interest to merit a 
more detailed investigation (including detailed rotation- 
curve fitting, lensing calculations, predictions concern- 
ing the formation of structure, cosmological predictions, 
etc.). We hope our summary of these constraints will fur- 
ther stimulate thinking about the relation between astro- 
physical puzzles and current ideas in high-energy physics. 
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